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Abstract 

During its first 10 years of orbital operations Swift dedicated approximately 11% of its observing time to blazars, carrying 
out more than 12,000 observations of ~ 1,600 different objects, for a total exposure time of over 25 million seconds. In this 
paper I briefly discuss the impact that Swift is having on blazar multi-frequency and time-domain astrophysics, as well as how 
it is contributing to the opening of the era of multi-messenger astronomy. Finally, I present some preliminary results from a 
systematic analysis of a very large number of Swift XRT observations of blazars. All the ’’science ready” data products that 
are being generated by this project will be publicly released. Specifically, deconvolved X-ray spectra and best fit spectral pa¬ 
rameters will be available through the ASDC ”SED builder” tool (https://tools.asdc.asi.it/SED) and by means of interactive tables 
(http://www.asdc.asi.it/xrtspectra). Innovative data visualisation methods (see e.g. http://youtu.be/nAZYcXcUGW8) are also being 
developed at ASDC to better exploit this remarkable and rapidly growing data set. 

Keywords: Active Galactic Nuclei; black hole physics; BE Lac objects; radiation mechanisms: non-thermal; 


1. Introduction 

Blazars are a rare type of AGN that are well known 
for their extreme observational properties such as super¬ 
luminal motion and highly variable non-thermal emis¬ 
sion across the entire electromagnetic spectrum, from 
radio waves to the highest energy y-rays. Their unique 
properties are attributed to emission of radiation by en¬ 
ergetic particles that move towards us in a magnetic field 
within a relativistic jet that happens to be pointing close 
to our line of sight d [2l . 

Blazars come in two flavours: Flat Spectrum Radio 
Quasars (FSRQs ) and BL Lacertae objects (BL Lacs). 
The main difference is in their optical spectra, with the 
former displaying broad emission lines and the latter 
being instead characterised by either a featureless con¬ 
tinuum, or by a spectrum that displays only absorption 
features (usually from the host galaxy) or weak narrow 
emission lines m 
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Figure 1. The SED of the blazar 3C454.3 assembled with over 60,000 
flux measurements obtained with many space and ground-based ob¬ 
servatories. The contribution of all Swift instruments (UVOT, XRT 
and BAT) is highlighted by the three grey vertical bands. 
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The radio to y-ray Spectral Energy Distribution 
(SED) of blazars typically displays two broad humps 
when plotted in y/(y) vs y space (see e.g. fig. 

It is widely accepted that the radiation associated to 
the low-frequency hump in the SED of blazars is due to 
synchrotron emission from relativistic electrons mov¬ 
ing in a magnetic field within the jet. The nature of 
the high energy emission, may instead be attributed to 
two intrinsically different mechanisms, described within 
scenarios generally referred to as leptonic and hadronic 
models n In leptonic scenarios the y-ray emission is 
assumed to be due to inverse Compton radiation EllSI 
whereas in leptohadronic models it could be due to pro¬ 
ton synchrotron radiation cmi or it may have a photo- 
hadronic origin a 

The peak energy of the synchrotron component 
(Vpeak) ranges from about 10^^-^ Hz (equivalent to ap¬ 
proximately 0.01 eV) to over 10^^-^ Hz (~13 KeV), re¬ 
flecting the maximum energy at which particles can be 
accelerated 0. Sources where 
Hz (-0.4 eV) are called low synchrotron peaked blazars 
(LBL), while those where > 10^^ Hz (-4 eV) are 
called high energy peaked blazars, or HBL ifTOl . 

Eor a long time blazar research has been confined to 
a small community of specialists. This is rapidly chang¬ 
ing as these sources are receiving much more general 
attention after they have been found to be the largest 
population of extragalactic sources in microwave (be¬ 
tween -30 and - 200 GHz), and in y-ray (0.1< E < 
100 GeV) surveys, as well as the most common extra- 
galactic objects appearing in catalogs of TeV detected 
sources. Most of the blazars found in microwave sur¬ 
veys are of the LBL type (both ESRQs and BE Lacs), 
while HBL BE Lacs dominate the y-ray and the VHE 
(E >100 GeV) extragalactic sky. 

Blazars are powerful extragalactic sources capable of 
accelerating particles to very high energies and there¬ 
fore they are considered as prime candidates for multi¬ 
messenger astrophysics. 

2. Swift and blazars 

During its first ten years in orbit the Swift obser¬ 
vatory mi carried out more than 12,000 observations 
of approximately 1,600 blazars, for a total of 25.3 
Ms of X-Ray Telescope (XRT, lIT^ ) net exposure 
time, or approximately 11% of the entire scientific 
program. About 50% of the objects listed in the most 
up to date catalog of blazars (BZCAT, 5* edition 
ini) have been observed at least once. Some of 
the well known and best studied objects have been 



Time [years] 

Figure 2. The rate of Swift observations including at least one blazar 
in the XRT field of view. Note the rapid increase from the beginning 
of the scientific operations in early 2005 up to 2010, after which date 
the rate decreased somewhat stabilising around an average level of 
^1500 observations/year. 


observed hundreds of times. Examples are the BE 
Lacs MRK 421 with 635 observations, MRK501 with 
412 observations, and the ESRQ 3C454.3 with 419 
pointings. Eig. illustrates how the rate of Swift 
blazar observations evolved in time, showing a rapid 
increase between 2005 and 2009, followed by an 
approximately constant level of 1,500 observations per 
year. This is quite a large number of pointings reflecting 
a scientific strategy which reserves to blazars a size¬ 
able fraction of the Swift overall observing program. 
The full data set is available from the official Swift 
archives in the U.S. (http://swift.gsfc.nasa.gov/), 
Italy (http://swift.asdc.asi.it/), and U.K. 
(http://www.swift.ac.uk/). However, current archive 
facilities only provide basic data sets and software 
packages for the data reduction. Scientific analysis 
must be carried out by the user by means of mission 
specific and additional analysis software. This is usu¬ 
ally done for a single exposure or for a limited number 
of observations; clearly normal users could hardly 
manage to analyse the entire data set. To facilitate 
the use of this valuable public archive at ASDC, also 
as part of the Swift team, we started a program to 
systematically analyse Swift blazar data to provide 
higher level ’’science ready” products and results, such 
as deconvolved spectra and best fit parameters for 
different spectral models. 

An example of the remarkable contribution that Swift 
is providing to multi-frequency studies of blazars is 
given by the SED of 3C454.3 shown in Eig. where 
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Figure 3. The first few lines of the interactive table through which the results of the ASDC XRT spectral analysis of blazars will be available on-line 
at http://www.asdc.asi.it/xrtspectra/ 


the simultaneous data obtained from the three Swift in¬ 
struments, UVOT, XRT and BAT, during the 419 obser¬ 
vations carried out before the end of 2014, is highlighted 
by three vertical bands. 

3. Swift XRT Data Analysis 

At the ASI Science Data Center (ASDC) we have 
processed in a uniform way, using the latest reduction 
software and calibration available, all XRT observations 
of blazars pointed by Swift at least 10 times (approxi¬ 
mately 120 blazars). The processing was carried out 
both on the overall data collected during each observa¬ 
tion and also on a orbit by orbit basis, thus allowing a 
time resolution of the order of 1.5 hours. 

The results of the analysis of approximately 5,000 
XRT observations (and nearly 9,000 orbits) will be pub¬ 
lished shortly Ga. The high level data products, such as 
deconvolved spectral data, light-curves, best fit parame¬ 
ters etc, will be available on-line within the ASDC SED 
builder tool (https://tools.asdc.asi.it/SED) and through 
interactive tables (http://www.asdc.asi.it/xrtspectra), a 
preview of which is shown in fig. Best fit spectral pa¬ 
rameters are given for the most commonly used spectral 
models (power law and log parabola), for each observa¬ 
tion and for every orbit during which sufficient statistics 
is collected for a spectral analysis. UVOT and BAT data 
for the same blazar sample will also be processed and 
the results will be published in a similar way in the fu¬ 
ture. In the following we show some preliminary results 
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Figure 4. The distribution of spectral slopes of the blazars most ob¬ 
served by Swift. FSRQs are shown in the upper panel, BL Lacs in the 
lower panel. 
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Figure 5. The correlation between X-ray spectral slope and X-ray flux 
in the BL Lac object MKN421. Note the very large range in flux (a 
factor of ^ 500) and spectral slope. More than 600 Swift observations 
are shown in this plot. 

of this work, which will lead to the establishment of one 
of the largest existing databases of high level scientific 
products and results on blazars. 

One example of these results is illustrated in fig 
where the distributions of the best fit power law spec¬ 
tral indices are shown for the subsamples of FSRQs and 
BL Lacs. Fig. shows instead how the spectral slope is 
related to source intensity for the case of MRK 421. The 
tight correlation apparent from this plot is well known 
[e.g. [151 US! 5 however the very large amount of data 
plotted in fig. [^provides unprecedented statistics both 
in spectral slope and in fiux range, this last spanning a 
factor of ^^500, probably a record value for soft X-ray 
variability of AGN. 

4. Time domain data 

The remarkable number of measurements accumu¬ 
lated over the last several years, thanks to Swift and 
many other space and ground-based observatories, are 
enabling new approaches to the analysis of luminosity 
variability in all parts of the electromagnetic spectrum. 
In this contribution blazars are used to provide some ex¬ 
amples of how this new opportunity is being exploited. 

Fig. shows the light-curve of 3C454.3 in the y- 
ray. X-ray, optical and radio band. It has been built us¬ 
ing the same data set used to assemble the SED of fig. 
[2 A simple visual inspection of this plot reveals that 
the light-curves at different frequencies share a simi¬ 
lar overall behaviour, but also show large differences. 
One technique oftenused to quantify these differences 



Figure 6. Multi-frequency light-curve of the blazar 3C454.3. In this 
plot the same y-ray, X-ray, optical and radio data extracted from the 
data set used for the SED of fi£] is used to illustrate the different 
behaviour in this source at different energies. 

is the discrete cross-correlation function (DCCF) Km 
between the light-curves at different energies. Fig. [7] 
shows the DCCF of 3C454.3 obtained by comparing 
the y-ray light-curve (1 GeV from Fermi-LAT observa¬ 
tions) with data taken in different energy bands, namely 
1 keV, 1 mm, 37 GHz and 15 GHz. As can easily be seen 
in fig. [7] the amount of correlation with y-rays is max¬ 
imum for 1 mm data, with no significant time lag. All 
other light-curves compared to IGeV data show a lower 
lever of correlation and significant time lags, ranging 
from about one month at IKeV to several months in the 
radio band, demonstrating the complexity of the emis¬ 
sion processes in blazars, which are subject to different 
dynamical timescales in various parts of the electromag¬ 
netic spectrum. This implies that while the use of si¬ 
multaneous data, often obtained through complex obser¬ 
vational campaigns involving manyground and space- 
based observatories (see e.g. CSl), is desirable and 
probably necessary, it is not sufficient for a full under¬ 
standing of the physics of blazars. 

A completely new technique to visualise variability 
in the SED of cosmic sources is to run in a sequence 
a set of frames, each representing the SED in a par¬ 
ticular time interval, just like in a movie. Of course, 
this can be done only if a sufficiently large number 
of multi-frequency measurements are available at all 
times. This is a requirement that is quite demanding 
but that has been already met by some of the bright 
and well known sources. As an example, fig. shows 
two frames of such a movie, which was built combining 
over 60,000 multi-frequency fiux measurements of the 
blazar 3C454.3. The full movie can be accessed on-line 
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Figure 7. The Discrete Cross Correlation Function on multi-frequency 
data of 3C345.3 illustrating how different energy bands (15 GHZ, 37 
GHz, 1mm, and 1 keV) show different time delay with respect to the 
y-ray emission at 1 GeV. 


at http://youtu .be/nAZYcXcUGW8 . 


5. Blazars and multi-messenger astrophysics 

Due to their large luminosity and SEDs often reach¬ 
ing the highest observed y-ray energies, blazars are the 
most powerful and most energetic cosmic accelerators 
known. 

An estimate of the maximum energy at which par¬ 
ticles are accelerated within this scenario can be ob¬ 
tained by measuring the peak (in y/(y) space) of the 
synchrotron hump, y^^^. Swift is sensitive in the soft 
and hard X-ray bands, where the of the most ener¬ 
getic blazars is located. Having observed about 50% of 
the known blazars the Swift archive provides the largest 
potential for the measurement of the synchrotron peak, 
and therefore the maximum energy at which particles 
can be accelerated in blazars. 

5. 1 . Neutrinos from blazars ? 

Since their detection at very high y-ray energies by 
Cherenkov telescopes, HBL blazars (which are the type 
of blazar most frequently detected in the VHE band) 
have been considered as candidate sources of extra- 
galactic neutrinos |[T9l[20l . 

The accuracy of the arrival direction of neutrinos as 
measured by current detectors, however, can be as large 
as several degrees, so a pure positional matching is 
hardly sufficient to guarantee a statistically sound asso¬ 
ciation. Padovani & Resconi (2014) ED suggested that 
seven BE Lac objects of the HBL type could be related 




Log(v) [Hz] 


Figure 8. Two frames extracted from the SED movie of the blazar 
3C454.3 that can be watched at http://youtu.be/nAZYcXcUGW8. 
The red points are simultaneous multi-frequency measurements taken 
within the time interval marked by the red dashed lines in the top part 
of each frame showing the IGeV light-curve. Green coloured text 
represents ground-based telescopes, orange text indicates satellite ob¬ 
servatories. The top frame shows the source in a state of growing 
intensity in which the emission is very large at all frequencies, with 
the exception of the radio and microwave bands. The lower frame il¬ 
lustrates an intensity state where the source is in a decreasing phase 
of its overall activity, but it is still bright at radio and microwave fre¬ 
quencies. 
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Figure 9. The SED of MKN 421 plotted together with the flux ex¬ 
pected from the twin-y-ray photons associated to IceCube neutrino 
with id 9, assuming that the association reported in ED is correct. 


to as many neutrino events recently reported by the Ice- 
Cube collaboration l22l on the basis of both positional 
and energetic diagnostics. That is, not only candidate 
blazar counterparts are required to be within the posi¬ 
tional uncertainly region of the IceCube neutrinos, but 
also their y-ray intensity, which is assumed to be tightly 
connected to neutrino emission, must be comparable to 
the observed neutrino flux. 

Fig. shows the SED of the BL Lac object MRK 
421, (one of the sources reported in ||2T|) including the 
expected flux from the twin y-ray photons assuming that 
the IceCube neutrino with ID 9 is indeed associated to 
this bright BL Lac object. Petropoulou and collabora¬ 
tors 1^ recently showed that the SED of this object 
can be interpreted in terms of a one-zone lepto-hadronic 
model. Assuming that this model is correct Padovani 
et al. 2015 (241, have estimated the total neutrino cos¬ 
mic background from the entire population of blazars 
and compared it with IceCube data. The expected over¬ 
all extragalactic intensity is close to the observed one at 
energies above ~ IPeV, supporting the hypothesis that 
blazars as a class are a major producer of extragalactic 
neutrinos. 

Other authors (251 have also argued in favour of a 
possible association of radio-brightblazars (mostly of 
the LBL type) in the TANAMI sample with IceCube 
neutrinos. Also in this case this hypothesis is based on 
theoretical predictions rather than on a signiflcant statis¬ 
tical signal. 

Upcoming deeper IceCube observations are clearly 
necessary to confirm or disprove the blazar-neutrino as¬ 
sociation. 


5.2. Could blazars be the sources of UHECRs? 

The origin of Ultra High Energy Cosmic Rays (UHE¬ 
CRs) is a major question that after many years of obser¬ 
vations still remains unanswered (26l . Many such parti¬ 
cles with energy exceeding 50 EeV have been detected 
by the most sensitive cosmic ray observatories such as 
the Telescope Array in the northern hemisphere and the 
Pierre Auger observatory in the south. Their arriving di¬ 
rections are not correlated with the Galactic coordinates 
and their energy is so large that they are not expected 
to be deflected much by magnetic fields therefore re¬ 
taining the information about the source of origin. High 
expectations have therefore been expressed in the liter¬ 
ature for the association of UHECRs with extragalactic 
objects. Attempts have been made to associate AGN 
to UHECRs (27l[28l by cross-matching the positions of 
nearby AGN in the Veron-Cetty & Veron (2006) catalog 
(29l and of AGN detected in the Swift-BAT surveywith 
Auger UHECRs finding what it seemed to be a likely 
association supported by an encouragingly high confi¬ 
dence level. The search was limited to sources closer 
than -100 Mpc due to the fact that extragalactic UHE¬ 
CRs are expected to be severely attenuated by pion pho¬ 
toproduction interactions with the cosmic microwave 
background (CMB) radiation, the so called GZK effect. 
This high level of statistical confidence, however, was 
not confirmed by subsequent searches carried out when 
larger samples of UHECRs became available (30l[3T]| . 
The disappointing outcome curtailed the high expec¬ 
tations for the discovery of the astrophysical counter¬ 
parts of UHECRs. In fact, these initial expectations that 
might have been resting on optimistic confidence lev¬ 
els possibly estimated without fully taking into account 
of the penalty associated with the many attempts per¬ 
formed with different cuts and on different samples in 
the search for the long-awaited correlation with astro- 
physical sources. 

These searches were however limited to mostly ra¬ 
dio quiet AGN, that is sources that do not radiate at 
very high energy. Blazars, particularly those with high 
energy synchrotron peak (HBL) that are often detected 
at VHE energies, are the extragalactic sources that are 
known to accelerate particles to the highest observed en¬ 
ergies, and therefore are natural candidates as sources of 
UHECRs. In this context Tinyakov & Tkachev (2001, 
El) searched for and found a high confidence correla¬ 
tion between BL Lacs and UHECRs. At ASDC in the 
months before the Swift 10 year meeting, we have tested 
this hypothesis by comparing the arriving directions of 
UHECRs from the TA and Auger samples with the po¬ 
sition of the IWHSP blazar sample (33l, the most com¬ 
plete and largest sample of blazars of the type that are 
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Figure 10. AitofF projection in Galactic coordinates of the UHECRs 
detected by the TA and Auger cosmic ray observatories, together with 
the positions of IWHSP blazars that could be detected by the cur¬ 
rent generation of Cherenkov telescopes. Large circles highlight the 
matches between UHECR and blazars. 


expected to emit at VHE energies. In particular we used 
the subset of 110 IWHSP blazars that are bright enough 
to be detected by the current generation of Cherenkov 
telescopes, and cross-matched it with the arrival direc¬ 
tions of TA and Auger UHECRs obtaining a combined 
probability of chance occurrence lower than 10“^ |[34l . 
However, a preliminary verification of this encouraging 
result on the enlarged sample of UHECRs presented by 
the Auger collaboration in late 2014 oa once again did 
not confirm the significance of the correlation. To il¬ 
lustrate the status of the search at the time of writing 
fig. plots the arrival directions of the UHECRs from 
TA (red points) and Auger (blue and green circles) in 
Galactic coordinates, together with the positions of the 
IWHSP sources that are expected to be above the sensi¬ 
tivity threshold of the current generation of TeV tele¬ 
scopes (grey points). The locations where there is a 
matching between blazars and UHECRs (within 3.2 de¬ 
grees for the case of TA and 2.5 degrees for Auger) are 
indicated by open large circles. 

Work is still in progress as we are currently using new 
Swift data to identify the sources with the largest syn¬ 
chrotron peak energy and therefore refine the sample by 
selecting the closest, most powerful and most energetic 
objects. The results will be presented in a future publi¬ 
cation. 


6. Conclusion 

Since the beginning of its orbital operations Swift 
played a major role in accumulating what has now 
become one of the largest public archives of multi¬ 
wavelength data covering a time span of over ten years 
and including more than 12,000 observations of 1,600 


blazars. This is certainly one of the richest existing 
multi-frequency open-access archive that is easily ac¬ 
cessible through on-line services. 

The large potential for new discoveries and what can 
be done with this powerful archive has been illustrated 
with a number of examples. 

In this paper I have also shown how, thanks to Swift 
and to many other space and ground-based astronomical 
facilities, innovative techniques are emerging enabling 
significant progress in the analysis of multi-frequency, 
multi-temporal astronomical data. These new tech¬ 
niques exploit the rich content of the available datasets 
in a deeper and more complete way compared to pre¬ 
vious approaches, and will no doubt lead to discover¬ 
ies that would not possible through a conventional ap¬ 
proach, thus making archival research an increasingly 
effective tool for modem astrophysics. 

Multi-messenger astrophysics is the additional fron¬ 
tier. Although an irrefutable proof that blazars are con¬ 
nected to neutrinos or UHECRs has yet to be found, 
these peculiar objects firmly remain among the most 
promising extragalactic sources holding the potential 
for opening soon this new and exciting window on the 
Universe. 
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